KtKUK  1  UUUUMtN  1  A  1  IUIN  KAtit  |  OMB  Na  0704-0188 

Public  reporting  burden  for  this  ejection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  data  sources, 

gathering  and  njaintevvng  the  data  needed,  and  completing  and  reviewing  the  collection  of  information  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection 

of  information,  including  suggestions  for  reducing  this  burden  to  Washington  Headquarters  Service,  Directorate  for  Information  Operations  and  Reports, 

1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget, 

Paperwork  Reduction  Project  (0704-0188)  Washington,  DC  20503. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 

1 .  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE 

1  Jan  2005  Final  Technical  Report 

3.  DATES  COVERED  (From -To) 

1  Jan  05-31  Dec  08 

4.  TITLE  AND  SUBTITLE 

MODELING  A  HALL  THRUSTER  FROM  ANODE  TO  PLUME  FAR 
FIELD 

5a.  CONTRACT  NUMBER 

FA9550-05-1  -0042 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 
lain  D.  Boyd 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Department  of  Aerospace  Engineering 

University  of  Michigan 

Ann  Arbor,  MI  48109-2140 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRE 

SS(ES) 

USAF/AFRL 

AFOSR 

875  North  Randolph  Street 

Arlington  VA  22203 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 
AFOSR 

11.  SPONSORING/MONITORING 

AGENCY  REPORT  NUMBER 

rt^-ftteVft-Tg-ZPia-  0513 

12.  DISTRIBUTION  AVAILABILITY  STATEMENT 


Distribution  Statement  A:  Approved  for  public  release.  Distribution  is  unlimited. 


13.  SUPPLEMENTARY  NOTES 


14.  ABSTRACT 

Two  dimensional  ax  symmetric  simulations  of  xenon  plasma  plume  flow  fields  from  a  D55 
Anode  layer  Hall  thruster  is  performed.  A  hybrid  particle-fluid  method  is  used  for  the 
Simulations.  The  magnetic  field  surrounding  the  Hall  thruster  exit  is  included  in  the 
Calculation.  The  plasma  properties  obtained  from  the  hydrodynamic  model  are  used  as 
boundary  conditions  for  the  simulations.  In  this  simulation,  the  Boltzmann  model  and  a 
detailed  fluid  model  are  used  to  compute  the  electron  properties,  the  direct  simulation 
Monte  Carlo  method  models  the  collisions  of  heavy  particles,  and  the  Particle-ln-Cell 
method  models  the  transport  of  ions  in  an  electric  field.  The  accuracy  of  the  simulation  is 
assessed  through  comparison  with  various  measured  data.  It  is  found  that  a  magnetic  field 
Significantly  affects  the  profile  of  the  plasma  in  the  Detailed  model.  For  instance,  in  the  case 
of  zero  magnetic  field,  the  plasma  has  a  potential  about  80  V  at  10  mm  from  the  thruster 
exit,  while  in  the  case  of  a  magnetic  field  included,  the  plasma  potential  is  about  60  V. 
Results  predicted  by  the  Detailed  model  with  the  magnetic  field  are  found  to  be  in  better 
agreement  with  experimental  data. 

15  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

Unclassified 

18.  NUMBER 

OF  PAGES 

11 

19a.  NAME  OF  RESPONSIBLE  PERSON 

a.  REPORT 

Unclassified 

b.  ABSTRACT 

Unclassified 

c.  THIS  PAGE 

Unclassified 

19b.  TELEPONE  NUMBER  ( Include  area  code) 

Standard  Form  298  (Rev.  8-98) 

^^rescn&ed  oy  Ahtal-bld 

MODELING  A  HALL  THRUSTER 
FROM  ANODE  TO  PLUME  FAR  FIELD 


AFOSR  GRANT  FA95 50-05- 1-0042 


Iain  D.  Boyd 

Department  of  Aerospace  Engineering 
University  of  Michigan 
Ann  Arbor,  MI  48109-2140 


Abstract 


Hall  thrusters  represent  an  efficient  form  of  electric  propulsion  that  is  being  implemented 
on  USAF  satellites  for  station-keeping  tasks.  Numerical  simulation  of  the  plasma  flow 
inside  Hall  thrusters  is  playing  an  increasingly  significant  role  in  the  development  of  the 
thrusters  to  better  understand  propulsion  and  lifetime  performance.  Numerical  simulation 
of  the  plumes  of  Hall  thrusters  is  required  to  accurately  assess  spacecraft  integration 
issues.  The  work  being  pursued  here  seeks  to  develop  for  the  first  time  an  end-to-end 
Hall  thruster  simulation  capability  that  begins  with  propellant  injection  at  the  thruster 
anode,  and  ends  in  the  plume  far  field.  The  development  of  a  comprehensive  simulation 
capability  is  critical  for  a  number  of  reasons.  The  main  motivation  stems  from  the  need 
to  directly  couple  simulation  of  the  plasma  discharge  processes  inside  the  thruster  and  the 
transport  of  the  plasma  to  the  plume  far  field.  The  simulation  strategy  will  employ  two 
existing  codes,  one  for  the  Hall  thruster  device  and  one  for  the  plume.  The  coupling  will 
take  place  in  the  plume  near  field  region  that  has  not  previously  been  modeled  in  detail. 
The  simulation  process  will  be  assessed  through  application  to  Hall  thrusters  for  which 
measured  data  exist. 

Recent  accomplishments 

Significant  progress  has  been  made  in  improving  and  extending  the  modeling  capabilities 
for  both  Hall  thruster  devices  and  Hall  thruster  plumes.  In  particular,  we  have  focused  on 
thruster  and  plume  modeling  of  the  D55  anode  layer  Hall  thruster  due  to  the  wide 
availability  of  detailed  flow  field  and  performance  measurements.  In  the  following 
sections,  we  summarize  this  progress. 

Hall  Thruster  Channel  Modeling 

Hall  thruster  device  modeling  efforts  were  focused  on  calculating  the  thruster  exit 
condition.  A  model  for  the  thruster  channel  is  developed  [1],  This  model  is  used  to 
simulate  the  flow  within  the  channel  of  the  D55  TAL  thruster. 

We  consider  three  conditions  corresponding  to  three  different  experiments.  These 
experiments  conducted  at  the  University  of  Michigan  [2]  [3],  TsNIIMASH  [4],  and  the 
University  of  Tennessee  Space  Institute  (UTSI)  and  Lockheed  Martin  Astronautics 
(LMA)  [5],  It  is  reported  [6]  that  some  portion  of  the  plasma  plume  of  a  D55  thruster 
consists  of  doubly  charged  ions.  In  this  research,  the  number  fraction  of  double  xenon  ions 
is  assumed  to  be  0.2. 

Contour  plots  of  the  plasma  density  and  potential  fields  are  shown  in  Figs.  1  and  2, 
respectively.  The  anode  is  on  the  left  side  of  the  figure,  while  the  exit  plane  is  at  the  right 


border.  The  white  spaces  above  and  below  the  channel  in  the  figures  represent  the  outer 
and  inner  wall  thicknesses,  respectively,  in  relation  to  the  channel  domain.  The  velocity 
distributions  are  based  on  the  ion  temperature,  which  can  be  assumed  to  be  the  same  as 
the  initial  neutral  temperature,  which  in  turn  is  based  on  a  wall  temperature  of  1000  K. 
The  velocity  distributions  are  used  to  create  a  distribution  of  ion  flow  angles. 

The  radial  distribution  of  the  ion  number  density  and  velocity  at  the  thruster  channel  exit 
are  used  for  the  plasma  plume  simulation. 


Figure  1 :  The  simulated  plasma  density  Figure  2:  The  simulated  plasma  potential 
field  within  the  D55  TAL  thruster  field  within  the  D55  TAL  thruster 


Hall  Thruster  Plume  Field  Modeling 

Our  recent  studies  have  focused  on  careful  assessment  of  our  existing  Hall  thruster  plume 
models  using  the  plume  data  measured  for  a  D55  TAL  thruster  [2][3][4][5]. 
Experimental  investigation  of  the  magnetic  field  distribution  near  a  Hall  thruster  shows 
that  the  radial  magnetic  field  component  is  much  larger  than  axial  component  [2].  Thus  for 
simplicity,  only  the  radial  component  of  the  magnetic  field  is  considered. 

Plasma  Potential 

Figures  3a  and  3b  show  radial  profiles  of  plasma  potential  at  axial  distances  of  10mm  and 
500mm  from  the  thruster  exit  plane,  respectively.  At  10mm,  the  Boltzmann  model 
overpredicts  the  potential.  The  Detailed  model  captures  the  shape  quite  well  although  it 
overpredicts  the  potential  too.  However,  the  result  is  better  than  the  Boltzmann  model 
and  if  the  magnetic  field  is  considered  the  improvement  is  much  better.  At  a  distance  of 
500mm  from  the  thruster,  the  Boltzmann  model  again  greatly  overpredicts  the  potential 
whereas  the  Detailed  model  reproduces  fairly  well  the  measured  profiles. 

Ion  Current  Density 

Ion  current  density  profiles  predicted  by  the  simulation  are  compared  with  the 
experimental  data  in  Figs.  5a  and  5b  along  radial  lines  located  at  10  and  500mm  from  the 
thruster  exit  plane,  respectively.  Both  the  Boltzmann  model  and  the  Detailed  model  give 


good  prediction  at  10mm.  At  500mm,  the  Detailed  model  shows  better  agreement  with 
measurements  though  both  models  underpredict  the  measured  values.  This 
underprediction  of  ion  current  density  over  the  entire  domain  suggests  that  the  simulation 
may  over  accelerate  ionized  particles  in  the  radial  direction.  This  feature  is  consistent  with 
the  comparisons  of  electron  number  density  shown  in  Figs.  6. 


a:  10  mm  from  the  thruster  exit  plane 


b:  500  mm  from  the  thruster  exit  plane 


Figure  4:  Radial  profiles  of  plasma  potential 
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Figure  5:  Radial  profiles  of  ion  current  density 

Electron  Number  Density 

Measurements  of  electron  number  density  are  compared  with  the  simulations  for  radial 
profiles  at  10mm  and  500  mm  in  Figs.  6a  and  6b,  respectively.  The  simulation  values 


represent  the  total  charge  density  obtained  from  the  number  densities  of  the  Xe+  and  Xe++ 
ions.  At  10mm,  most  of  the  simulation  data  underpredict  the  measured  values.  The  peak 
electron  number  density  measured  at  both  stations  is  more  than  double  the  total  charge 
density  assumed  in  the  simulations  at  the  thruster  exit  plane.  One  interpretation  of  these 
comparisons  is  that  the  axial  component  of  electric  fields  in  the  simulation  is  so  strong 
that  the  acceleration  of  ions  is  overestimated  in  the  axial  direction.  One  possible  way  to 
address  such  differences  between  the  model  and  the  measured  data  would  be  axial 
confinement  of  electrons  caused  by  the  magnetic  field.  It  was  shown  that  the  axial 
component  of  the  magnetic  field  is  much  smaller  than  the  radial  component  in  SPT-100 
thrusters  [7],  although  there  is  no  direct  evidence  that  it  is  true  in  D55  thrusters. 
Measurements  of  the  axial  component  of  the  magnetic  field  are  required  to  help  resolve 
this  issue,  and  no  such  data  exist  as  of  now.  At  500mm,  the  Detailed  model  shows  better 
agreement  with  the  measurements  though  both  models  still  underpredict  the  measured 
values  over  the  entire  radial  profile. 


a:  10  mm  from  the  thruster  exit  plane  b:  500  mm  from  the  thruster  exit  plane 

Figure  6:  Radial  profiles  of  electron  number  density 

Electron  Temperature 

Figures  7a  and  7b  show  radial  profiles  of  electron  temperature  at  distances  of  10mm  and 
500mm  from  the  thruster,  respectively.  Electron  temperature  is  constant  in  the 
Boltzmann  model  and  here  we  show  the  value  of  lOeV  which  is  used  in  the  simulation.  At 
10  mm,  it  is  clear  from  the  measurement  that  there  is  significant  spatial  variation  in  the 
electron  temperature  caused  by  the  dynamics  of  the  plasma  as  it  enters  the  acceleration 
channel.  In  general,  although  the  Detailed  model  provides  reasonable  agreement  with  the 
measurements,  the  radial  gradients  predicted  by  the  model  are  smaller  than  the  measured 
data  indicate.  These  disparities  between  the  models  and  the  measurement  indicate  that 
more  elaborate  thruster  exit  boundary  conditions  are  needed.  In  fact,  it  is  known  that  the 
Detailed  model  is  relatively  more  sensitive  to  boundary  conditions  than  the  Boltzmann 
model.  Fig.  7a  shows  that  the  electron  temperature  decrease  if  the  magnetic  field  is 


considered.  If  the  magnetic  field  is  considered,  the  resulting  thermal  conductivity 
coefficient  is  decreasing  [8].  Therefore,  the  electrons’  motion  is  confined  and  the  electron 
temperature  decreases.  Fig.  7a  shows  that  the  Detailed  model  with  a  magnetic  field 
underpredicts  electron  temperature.  This  suggests  that  we  need  more  accurate  physical 
models.  Fig.  7b.  clearly  shows  that  the  Detailed  model  gives  good  agreement  with  the 
measured  data  in  the  far  field. 


a:  10  mm  from  the  thruster  exit  plane  b:  500  mm  from  the  thruster  exit  plane 


Figure  7:  Radial  profiles  of  electron  temperature 
Axial  Velocity  of  Xe+ 

Figure  8  shows  the  axial  velocity  profiles  at  a  radial  position  of  27.5  mm  which  is  along 
the  thruster  channel  center. 


Figure  8:  Axial  components  of  ion  velocity  at  a  radial  position  of  27.5mm 


It  is  clear  that  the  Boltzmann  model  fails  to  produce  sufficient  ion  acceleration  in  the  near 
field  of  the  plume.  This  is  an  expected  result  because  the  plasma  potential  gradient  shown 
in  Fig.  3  is  not  enough  to  give  ions  significant  acceleration.  The  Detailed  model  predicts 
strong  ion  acceleration  in  the  near  field  region  and  rapidly  accelerates  the  ions  from  the 
thruster  exit  velocity  of  1 5  km/s  to  a  value  of  about  1 8  km/s  that  corresponds  to  the 
measured  data  and  the  results  becomes  slightly  better  when  the  magnetic  field  is 
considered.  The  simulation  result,  however,  overestimates  the  axial  velocity  at  z=l  mm. 
One  possible  explanation  is  that  there  would  be  a  lot  of  CEX  collisions  in  front  of  the 
thruster.  Because  of  CEX  collision,  slow  ions  are  created,  and  these  slow  ions  would 
make  lower  the  central  value  of  the  ion  velocity  distribution  function. 

Relevance/Transitions 

The  plasma  plume  code  MONACO-PIC  used  in  this  research  is  distributed  to  Ken  Tatum 
of  AEDC  where  it  is  used  to  model  Hall  thruster  plume  expansions  created  in  the  1 2  V 
vacuum  facility  at  AEDC. 
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